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The potential concept that is successful in classical electrodynamics should also be applicable to 
the nonlinear electromagnetic forces acting on matter. The obvious method of determining these 
potentials should be provided by Helmholtz's theorem. It is found, however, that the theorem fails in 
most practical instances. Other methods to find the potentials - as pursued in plasma physics - are 
examined and found to yield functions which depend on the chosen coordinate system. Thus they 
cannot be considered as invariant potentials from which physical forces may be derived. Practical 
consequences of these mathematical findings are discussed. 
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1 Introduction 

In classical electrodynamics the force density acting on a volume element, which carries charge 
and current, is derived from the divergence of Maxwell's stress tensor. We have the well known result 
[1] which may also be derived directly from the Lorentz force: 



Both the electric and the magnetic force are nonlinear, since the electromagnetic field itself depends 
on the sources p and j according to Maxwell's linear equations (in vacuo) : 
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It is convenient to introduce a scalar and a vector potential 




B = rotA 
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from which the fields may be derived. The ansatz © satisfies equations ^ and (O automatically. 
Upon substitution into ([2]) and ([4]) it leads to two second order equations which may be solved by 
known mathematical methods. 

In view of the success of the potential method it is somewhat surprising that textbooks do not 
make use of it by applying it directly to ((!]). This should be possible, since, according to Helmholtz's 
theorem [5], any vector field may be decomposed to the form: 

F = Vp + rotu, divu — (7) 

This way the force F is expressed in terms of the scalar potential p and the vector potential m. The 
advantage of the decomposition becomes obvious when we want to balance the electromagnetic forces 
with mechanical forces inside conductors. In plasma confinement physics, e.g., one would like to 
balance the mechanical pressure gradient with the magnetic force: 

J X B = Vp (8) 

This requires finding a magnetic field configuration such that the rotational part of the cross product 
vanishes. The electric force in ([I} is of less practical importance, since the charge distribution is 
usually not even known. Instead one calculates the electric field from Ohm's law and determines 
subsequently the charge density from (2). In this case, however, Ohm's law for a moving conductor 
in the form 

E + vxB^rj] (9) 

requires us also to balance the gradient part of the cross product with the gradient part of the electric 
field in steady state. 

From these examples it should be clear that it is worthwhile to explore the possibility of expressing 
a cross product as the sum of a gradient and a curl. In Section 2 we take a closer look at Helmholtz's 
theorem and find certain limitations of its validity. It appears that it cannot always be used to solve 
first order partial differential equations of the type ([5]) or (jH]) in a similar way as the ansatz ^ satisfied 
the first order equations ^ and ([5]). In Section 3 we examine the usual methods of solving first order 
equations of type ([8]) or ([9]). The surprising result is that the scalar potentials turn out to depend on 
the chosen coordinate system and can, therefore, not represent physical quantities such as pressure 
or electrostatic potential. In Section 4 we discuss some physical consequences of our findings. 



2 Restrictions on the validity of Helmholtz's theorem 

Helmholtz's theorem is proven 2J in the infinite domain by referring to the uniqueness theorem 
which applies to the solution of Poisson's equation. When we take the divergence of the defining 
equation ([7]) 

Ap^divF (10) 
we obtain a Poisson equation which has the unique solution 

1 [[[ div'F{x') ^3 , 

subject to the boundary condition p (oo) = 0. Similarly, taking the rotation of ([7]) one obtains the 
vector Poisson equation 

Au^-rotF (12) 



with the solution 



1 fffrot'F(x') , 



An 



subject to the boundary condition u{oo) = 6. Substitution of the potentials pT|) and into ([T]) 
should yield the force field F (x) again. 

At first sight this conjecture is not expected to be true for a simple reason: The solutions and 
p3)l are completely analogous to the solutions for the electromagnetic potentials which arise from 
Poisson equations when the ansatz ([6]) is substituted into ^ and (j4|) in the static case. Whereas 



the sources p and j are localized, the potentials and the fields derived from them are finite also in 
the region outside the sources. The force ((1]), however, is proportional to the sources and, hence, 
vanishing outside the conductors. It would be surprising, if integrals of the type ([TT|) and (fT^ would 
yield finite electromagnetic fields outside the source region, but vanishing electromagnetic forces in 
agreement with ([T]). In fact, this is not the case in general, but holds only, when the force vanishes 
on the surface of the conductor where charge and current density are still finite. 
In order to demonstrate this we consider the magnetic force in its decomposition 

jxB = \/p + rotu (14) 

and choose a toroidal geometry as sketched in Figure [TJ The force in the direction of the minor radius 



has two contributions: 



Figure 1: Axisymmetric toroidal solenoid 



Fr = j^} - (15) 

Firstly, we consider a magnetic field which points only in toroidal direction and vanishes at the 
boundary r — r^,: 



B,n — C— 

^ 2R 



, Rq — r cos -d 



Insertion into the static equation (4) gives the current density components: 

C r 

^0 R 

so that the magnetic force has the components: 

(J2 J. (^2 _ ^2-) 



Fr = 



F^^O 



Now we must take the divergence of this force to obtain the Laplacian of the scalar potential: 



d Rr^ (rt-r^) 



^loRr dr \ 2i?2 
Substitution into (fTTj) yields the potential: 
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Taking the rotation of (|18p one obtains with (|13p for the toroidal component of the vector potential: 
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Rather than trying to evaluate these integrals analytically, we have calculated them numerically and 
compared the result for the force component in the midplane Z — 

dp du^ 



Fr = 

with the analytic expression resulting from ^TE\i : 



OR dZ 



Fr 
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(22) 
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In Figure O the individual components of are depicted together with their superposition and the 
analytic expression ((23l) . The agreement between the numerical integrals ((22l) and the exact expression 
((23)) is perfect within numerical accuracy. It is also observed that the individual contributions in ([22]) 
are finite outside the torus (which extends from R = .5 to 1.5), but cancel perfectly when they are 
superimposed. This result confirms the validity of Helmholtz's theorem for the chosen case. 
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Figure 2: Radial magnetic force Fr = represented by potentials 



Next we consider the second term in ()15p which we choose finite on the torus surface. Assuming 
a magnetic field of the form 

Bo^^ , Br = (24) 
we obtain with (|4]) the toroidal component of the current density: 

_ C {Rq + R) 



and with (jl5p the magnetic force component: 



Fr = - 



C r {Rq + R) 



(26) 



Taking divergence and rotation of this expression we obtain with ([TT|) and ([T5)l the integral represen- 
tation of the potentials: 
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The radial force in the midplane may be derived from these integrals with and compared with 
the analytic expression: 



MO 



(29) 



In Figure [3] we show the result. Obviously, the force as derived from the potentials ((27|) and ([28] 
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Figure 3: Radial magnetic force Fr = — B,? represented by potentials 

bears little resemblance with the given force (note that the scale of the F^ - plot has been reduced 
by a factor of 2 in Fig. 3). Furthermore, the superposition of the derivatives of the potentials does 
by no means vanish outside the torus. As already mentioned, this behaviour is to be expected 
for integrals like (P7|) and (pS)) . in particular, when the force is finite on the boundary as in (|^ . 
This holds also for the electric force - the first term in ([1]) - when the body carries a total charge 
resulting in a finite boundary field. Helmholtz's theorem is apparently not applicable when the vector 



field is discontinuous which is a common property of the electromagnetic forces except in specially 
constructed cases like ([16]). 

Our investigation so far reveals that the Helmholtz formalism for determining the potentials from 
which the electromagnetic forces could be derived will fail in most practical instances. On the other 
hand, physical differential equations like ([5]) or ([§]) can only have a solution, when the electric potential 
or the pressure do exist. It is therefore of great importance to explore different methods which seek 
a solution for the potentials. The next Section will deal with standard methods being used in plasma 
physics for obtaining solutions to the potential problem of the electromagnetic forces. 



3 Standard methods for solving first order vector field equations in 
plasma physics 

3.1 Ohm's law in a plasma 

In an early paper of 1958 Kruskal and Kulsrud [Sj examined the conditions under which a plasma 
can be confined by a magnetic field. They realized that Ohm's law ([9|) leads in conjunction with 
Faraday's law and the assumption of a steady state to a special first order equation —B ■ = r/j-B, 
since the electric field should be expressible as the gradient of a scalar, if rot E = 0. For this type of 
equation they coined the term 'magnetic differential equation'. Newcomb [3] formulated a solvability 

criterion ^ j ■ B j B^ dl — 0, where the integration has to be carried out along closed magnetic 

field lines. Pfirsch and Schliiter made an attempt to solve the magnetic differential equation 
without taking reference to Helmholtz 's theorem, although in effect they were trying to determine 
the gradient part of ([9]). They used a magnetic model field in the geometry of Figure 1 where the 
poloidal field lines closed themselves on circles. This field satisfied Maxwell's equation ([5]), but not 
exactly Q and the force equilibrium condition ([8]). Nevertheless, they obtained an estimate for the 
plasma diffusion coefficient in a toroidal configuration. Later on Maschke [B] generalized their results 
for a true equilibrium configuration that was supposed to satisfy both ([4]) and ([8]). He was able 
to confirm the order of magnitude of the predicted plasma losses. The present state of the art to 
calculate collisional transport in tokamaks including the effect of trapped particles can be found, e.g., 
in H. 

Following the method originally adopted in [B] we assume that an axisymmetric toroidal field 
configuration exists (e.g. 'tokamak' [5]), where the poloidal magnetic field components may be derived 
from the toroidal component of a vector potential: 



The poloidal field lines follow the contours of the surfaces "0 = const which are nested, but not 
necessarily concentric circles as in Figure [TJ The inner product of with the magnetic field yields 
the magnetic differential equation: 

BrEr + BzEz^VJ-B-E^B^ = s (31) 

In steady state the poloidal electric field components may be expressed as the gradient of a potential, 
whereas the toroidal component E^ — U /R could be induced by a transformer that produces a loop 
voltage U. The authors of [31 IH [51 [5] assumed that ([3T|) has a solution which means that a potential 
exists such that one can write: 

This equation specifies the derivative of the potential along the poloidal field lines. Straightforward 
integration over the contours of constant flux yields the 'potential' 

^ = ^oW~J^^dl (33) 







where dl is a line element on a contour ip = const. Newcomb's criterion imposes the integrability 
condition 

i/'=const 

which must be satisfied by s in order to avoid a multi-valued potential. 

Although the described procedure yields a function /) in a 'flux coordinate system', one 
cannot be certain whether it represents a potential which exists independently of the chosen coordinate 
system. This would be guaranteed if we had obtained the potential from the solution of a Poisson 
equation like (jlip . but this was not possible, as we do not have - at this point - any information on 
the velocity field that determines the Laplacian of the potential according to ([9]): 

A(j)::^ B -rotv-v-rotB (35) 

There is only some knowledge on the current and field configuration which enters into s as defined by 
([31]) . Maybe this was the reason why the Helmholtz method was ignored in the analysis of plasma 
equilibria. 

In order to guarantee the existence of the potential in any coordinate system, we must satisfy the 
necessary and sufficient condition rot E = 0, or: 

OEr ^ dEz 
dZ dR 

We can apply this condition when we take the gradient of ([5T|) 



(36) 



OEr dEz dBR dBz ds 

(37) 

dER dEz dBR dBz _ ds 
Br^ + Bz^ + Er^ + Ez^ ^ 

and derive with (|36p and ([5]) two magnetic differential equations for the electric field components: 



The solutions are in analogy to 



Er = RBz [fR W + Ir (V, 0] , Ir = / 



_ f _d_ dl 



lo 

I 



Ez^RBr [fz W + Iz (V', /)] , ^2 = / ^ 1^;^ j (39) 



dR \Bz J IVV-I 



d ( s \ dl 



where Newcomb's condition must also be satisfied. 

Instead of cylindrical coordinates one could also have used spherical coordinates with R = 
r sm9 , Z = r cos 6', Br = B,. sin 6* + Bg cos 9, Bz = B^ cos 6 — Bgs'mO . Equation (|?T|) trans- 
forms then into: 

Br Er + BeEe^s (40) 



and the steady state condition rot E = becomes: 

dEr _ dirEe) 



(41) 



d9 dr 

Taking now the gradient of ([^Dl) one obtains with and (|¥T|l the magnetic differential equations for 
the field components in the form: 



with the solutions: 

I 

Er^RBe [fr m + Ir (V', 0] , Ir I) = 



1 d f sr\ dl 



rdr \Bg J |Vi/^| 

lo 

I 

Eo^RBr [fe W + le 0] , lo (Va = / ^ (43) 

lo 

The same expressions must result, if we transform the solutions p9p directly into spherical coordi- 
nates: 



or: 



Er sin e + Ee cos9 ^ R {Br cos 6* - sin 6* ) (-0) + Ir. (V", 0] 
Er COS 6* - £'e sin 6* =R {Br sin 6 + Bg cos 6) [fz (^A) + /z ("0, 0] 

Er=RBr sin COS + /fl + fz + I z] + R Be [cos" e{fz + Iz)- sin' {fn + Ir)] 



(44) 



(45) 



Ee^RBr [cos2 9 {fn + Ir) - sin' {fz + Iz)]-RBe sin 9 cos ^? [/« + Ir + fz + Iz] 
In order to make (|45p compatible with (|43)) it is apparently necessary to require: 

fR + lR + fz + lz = (46) 
On the other hand, if one substitutes the solution (|39p into (PT|) . one obtains: 

i? Sfl^ (/if. + /fl. + /z + /z ) - s (47) 



The inhomogeneous part s of the magnetic differential equation ([32)1 must apparently vanish, if the 
function is to exist as a potential independent of the coordinate system. This means, of course, 
that a scalar potential for the arbitrarily given vector field rjj — vxB in ^ does not exist. Only a 
function (j) = (pQ {ip) would satisfy Ohm's law (|32p when s — 0, but one could not consider it as a true 
potential either. This will be discussed in the following subsection. 

3.2 The force equilibrium in a plasma 

The stationary force balance in a magnetized plasma cannot be easily satisfied. In fact, there are 
not any analytic solutions of ([8]) known, except in the axisymmetric case |9l [10] . From the toroidal 
component of ([5]), namely jzBr — JrBz = 0, follows with ^ a homogeneous magnetic differential 
equation for the toroidal field component: 

B-V{RB^)^Q (48) 

and the inner product of ([5]) with the magnetic field leads to another homogeneous magnetic differ- 
ential equation 

B-Vp = (49) 

Inserting the obvious solutions: 

B^^F{^P)/R, p^p{^) (50) 

into the poloidal components of ([5]) one finds with ^ an expression for the toroidal component of 
the current density: 

'^^ dtp fio R dip 



Substitution into the toroidal component of ^ yields with ([50)1 the famous nonlinear Liist-Schliiter- 
Grad-Rubin-Shafranov equation [llj : 

It is the basis for calculating numerically axisymmetric plasma equilibria. Recently the equation has 
been modified [T^] to include approximately the effect of magnetic islands which break axisymmetry. 

One must expect, however, that a function p {ij)) which depends on a single vector component, 
namely Rcr x e<^, cannot represent a physical pressure that must be defined independent of the 
coordinate system. In order to see this we consider an analytical 'Soloviev solution' [9] 



^^t(coR' + bRl)Z^ + ^{R'-RlY , F^^Fl~2bRl^P (53) 

which satisfies (j52p by construction. The pressure is given as a linear function of if): 

p{R, Z)=pn- {a/iia) i'=po- (a/Aio) RA^ (i?, Z) (54) 

Transforming this equation into a Cartesian coordinate system yields: 

p{x, y, z) ^ Po - (a/no) {x Ay ~ y A^) (55) 

Rotation of the coordinate system around the y - axis by an angle a according to the transformation 
rules 



x' sin a + z' cos a , y — y' 



(56) 



Ax — Ax' cos a — Az' sin a , A^ = Ax' sin a + A^' cos a , Ay = Ay' 

results in a pressure field 

p{x', y', z) =po- {a/^J■o) [{x Ay' - y' Az')cosa+ {y A^' - z' ylj^- ) sin a] (57) 

which depends not only on the new coordinates, but also on the rotational angle a. It is, therefore, not 
a scalar in the usual definition: p' {x') — p{x) for x' ~ x. Our conclusion is then that the magnetic 
force j X B does not have a scalar potential which could be identified with a physical pressure. 
The so-called 'force-free' configuration 

jxB = 6 (58) 

which plays a role in very diluted astrophysical plasmas, deserves special mentioning. Obviously, it 
makes little sense to replace the - vector in ([55]) by a gradient and a curl. Whether a solution to 
(|58p exists, must be examined in a different way. In the axisymmetric case we have from (152^ : 



dF 

A*iP = -F—^-giiJ;) (59) 

Applying Stokes's theorem on this equation by integrating the toroidal current density over the area 
enclosed by a magnetic surface one has: 

IWI „ If , ,,dRdZ 
J rt #.9W^^ (60) 

With the inverse function 



R .If R 



one obtains from ([59)1 : 



V' = / ^ dF, W = - VF (61) 
J 9{F) 9 

/:,*F = -^--^\VF\'—(-\ (62) 

F F^ ^ dF\g) ^ ' 



Stokes's theorem applied on this equation gives: 



VF| „ Iff 9^ 5,^^,2 d fF\\dRdZ 



R Jj \F F' ' dF \ g J J R 

Substitution of ((6T|) into ((60)) yields on the other hand 



F f IVFI „ If dRdZ 

Elimination of the line integral over the poloidal current density on the left-hand-sides of ([55]) and 
(|64p . and using (j6ip again results in an integral equation 

which can only be satisfied for g = FdF/dip = 0. This may be demonstrated by choosing g = const 
and performing a partial integration on (j65p : 



Z(i^,^).i? r ^i^il)i^ + F#f^|^ + 4|V^r)^^ (66) 



i? F J R jf \ F^ dZ F^ ' ^' J R 

i/>=const i/'=const 

The first term on the right-hand-side cancels against the left-hand-side, whereas the second term 
cannot vanish on every surface except for g — 0. This may be verified by adopting the force- free 
Soloviev solution ([55]) with a = and g — — Wc must conclude then that equation (|58p in 
conjunction with Ampere's law Q does not have a solution. 



4 Discussion and conclusion 

Our attempts to determine scalar and vector potentials for the electromagnetic forces were only 
partially successful. Thanks to Helmholtz's theorem the task of formulating and calculating the 
potentials should be an easy exercise, but it was found that in most practical cases - where the forces 
are discontinuous across the boundary of the conductors - the theorem is not applicable. 

On the other hand, the potentials must exist, if the electromagnetic forces are to be balanced 
against mechanical forces in a stationary state, in particular in fluid conductors like plasmas, for 
example. We have examined the usual method for finding the electrostatic potential in Ohm's law 
and the pressure potential in the force balance of a plasma without having regard to Helmholtz's 
theorem. Although it appears possible to find solutions in specially chosen coordinate systems, the 
functions from which the forces may be derived cannot be identified with physical potentials, as they 
are not invariant against arbitrary transformations of the coordinate system. 

Our finding has important implications on the physics of magnetic plasma confinement. If it 
is not possible to find solutions for the steady state equations, a plasma can only be confined in 
a 'quasistationary' state at best. This means that the partial time derivatives neither in the law of 
induction, nor in the equations of motion can be neglected. The plasma may develop into a fluctuating 
turbulent state where strict axisymmetry does not prevail anymore. The observed 'anomaly' of plasma 
losses is probably due to this lack of a well balanced equilibrium in confining devices like tokamaks, 
stellarators, reversed field pinches etc. Harold Grad - one of the authors of equation ([52]) - came in 
his last paper [13] to a similar conclusion: "The discovered lack of pressure balance is not a special 
consequence of one idealized model but carries over broadly to a variety of physical refinements 
and generalizations. Indeed, the mathematical result has profound physical consequences which 
have gradually entered the field in the context of rational rotation number resonances and multiple 
helicities, island formation, and turbulence." 



References 

[1] J. D. Jackson, Classical Electrodynamics, Second Edition, (JohnWiley & Sons, Inc., New York, 
1975) Sect. 12, p. 607, eq. 12.141 

[2] G. Arfkcn, Mathematical Methods for Physicists, Third Edition, (Academic Press, Inc., Orlando, 
1985) Chapter 1.15 

[3] M. D. Kruskal, R. M. Kiilsrud, The Physics of Fhiids 1 (1958) 265 

[4] W. A. Newcomb, The Physics of Fluids 2 (1959) 362 

[5] D. Pfirsch, A. Schluter, Max-Planck-Institut fiir Physik und Astrophysik, Report MPI/Pa/7/62 

(1962) (unpubhshcd) 

[6] E. K. Maschke, Plasma Physics 13 (1971) 905 

[7] D. A. Gates, H. E. Mynick, R. B. White, Physics of Plasmas, 11 (2004) L45 
[8] John Wesson, Tokamaks, (Clarendon Press, Oxford, 1987) 

[9] S. Soloviev, Sov. Phys. JETP 26 (1968) 400; Zh. Eksp. Teor. Fiz. 53 (1967) 626 

[10] C. V. Atanasiu, S. Giinter, K. Lackner, I. G. Miron, Physics of Plasmas 11 (2004) 35 

[11] R. Liist, A. Schliiter, Zeitschrift fiir Naturforschung 12 (1957) 850 

V. D. Shafranov, Sov. Phys. JETP 6 (1958) 545; Zh.Eksp.Teor. Fiz. 33 (1957) 710 

H. Grad, H. Rubin, Proc. 2nd U. N. Int. Conf. on the Peaceful Uses of Atomic Energy Geneva 

1958, Vol. 31, 190, Columbia University Press, New York (1959) 

[12] X. Liu, J.D. Callen, C. G. Hegna, Physics of Plasmas 11 (2004) 4824L 

[13] H. Grad, International Journal of Fusion Energy 3 (1985) 33 



